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ABSTRACT 


The  objective  of  this  report  is  the  presentation  of  an  analytical- 
stochastic  model  capable  of  predicting  relevant  statistical  scattering 
features  of  electromagnetic  waves  propagating  within  vegetated  environ¬ 
ments.  The  propagation  phenomena  are  described  by  formulating  the 
scattering  associated  with  random  permittivity  fluctuations  superimposed 
on  a  lossy  deterministic  background  slab.  The  mean  backscattered  powej-, 
its  variance  and  one-point  distribution  are  calculated.  The  spectral 
characteristics  of  clutter  from  windblown  foliage  are  investigated  using 
two  models;  one  presuming  the  velocity  field  to  constitute  a  multi- variate 
normal  process  and  another  presuming  the  scatterer’s  motion  to  be  quasi¬ 
harmonic  over  limited  time  segments. 
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CLUTTER  RETURN  FROM  VEGETATED  AREAS 


INTRODUCTION 


The  objective  of  this  report  is  the  presentation  of  an 
analytical -stochastic  model  capable  of  predicting  relevant  statis¬ 
tical  scattering  features  of  EM  waves  propagating  within  vegetated 
environments.  Of  interest  are  such  features  as  the  spectral 
content  of  the  clutter  return  as  well  as  the  received  power  (and 
related  statistics)  from  gated  range-cells.  Other  pertinent  wave 
features  not  discussed  but  which  can  be  pursued  within  the  frame¬ 
work  of  the  suggested  model  are  mentioned  at  the  end  of  this 
section. 

The  forest  constitutes  a  complicated  transmission  mediunl 
Its  characterizing  parameters  (i.e.,  permittivity  and  conductivity) 
are  random  functions  of  space  and  time.  Furthermore,  its  inter¬ 
faces,  particularly  the  air-vegetation  interface  which  may  play  a 
major  role  in  the  transmission  process  (e.g.,  through  its  support 
of  a  "lateral"  wave)  are  "rough".  Further  discussion  of  the 
stochastic  character  of  the  interfaces  and  the  corresponding 
effects  on  the  EM  radiation  is  omitted  at  present.  This  by  no 
means  implies  lack  of  significance.  It  is  the  surface  roughne$s 
which  conceivably  causes  the  deterioration  of  the  "lateral"  wave, 
The  forest  constitutes  an  ensemble  of  various  scattering 
constituents.  Fortunately,  not  all  influence  the  EM  characteris¬ 
tics  in  the  same  manner.  The  outstanding  difference  between 


these  scatterers  is  their  physical  dimension  compared  to  the 
wavelength  of  the  incident  radiation. 

Small-scale  scatterers  are  characterized  by  small  scattering 
cross-section.  Their  relative  contributions  to  the  scattered 
field  (the  elementary  scattering  event  is  virtually  isotropic) 
are  expected  to  be  small.  The  presence  of  small-scale  scatterers, 
however,  is  extremely  significant  because  of  the  modifications 
they  impose  upon  the  effective  propagation  features  of  both  the 
incident  as  well  as  the  scattered  radiation.  Multiple-scattering 
events  experienced  by  either  the  incident  wave  or  by  the  wave 

(1-3') 

back-scattered  to  the  receiving  antenna  have  a  two-fold  effect v  ' 

(a)  Scattering  losses  (additional  to  the  ever-present 
absorption  losses)  reduce  the  efficiency  of  the  channel 
and  must  be  properly  accounted  for.  The  medium  which 
effectively  describes  the  propagation  of  the  mean  wave 
predicts  losses  which  may  become  considerably  higher  than 
those  anticipated  from  absorption  alone. 

(b)  In  addition  to  the  scattering  losses  the  mean  wave 
experiences  an  effective  slowdown.  The  medium  seems 
denser  than  that  expected  in  the  absence  of  scattering. 

The  small-scale  scattering  effects  indicated  above  are 

accounted  for  by  a  proper  choice  of  the  effective  background 
permittivity  and  conductivity.  The  relative  contributions  of 
small-scale  scatterers  to  the  scattered -wave  may  not  be  small 
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throughout  the  contributing  spectrum.  Specifically,  attention 
must  be  paid  to  the  spectral  "tail"  region  to  which  small-scale 
scatterers  may,  despite  their  small  cross  section,  contribute 
significantly  since  they  may  attain  relatively  high  velocities. 

The  effects  of  scattering  centers  contributing  to  the  bulk 
of  the  back-scattered  energy  cannot  be  accounted  for  in  an  analo¬ 
gously  simple  manner  since  the  scattering  is  controlled  by  the 
random  deviations  of  the  (complex)  permittivity  function  from  the 
background,  a  more  accurate  description  of  the  space-time  medium 
fluctuations  is  called  for.  First  order  scattering  theories  high¬ 
light  the  need  for  more  detailed  knowledge;  that  contained  in  the 
space-time  auto  and  cross  correlation  functions  of  the  pertinent 
electromagnetic  parameters. 

At  relatively  low  frequencies  (say  below  100  MHz)  virtually 
all  the  forest  constituents  are  in  the  small-scale  category.  The 
effect  of  the  scattering  on  the  propagation  characteristics  of  the 
mean  wave  is  totally  describable  by  an  effective  background  selec¬ 
tion  leading  to  the  familiar  deterministic  "forest-slab"  model 
The  scattered  wave  is  both  small  in  magnitude  and  inconsequential 
to  applications  at  these  frequencies  and  thus  has  been  totally 
ignored.  As  the  frequency  rises  to  the  range  which  includes  UHF 
and  L-band,  the  scattered  wave  can  no  longer  be  ignored.  The 
small-scale  scattering  constituents  cannot  be  ignored  either; 
their  determination  of  the  background  slab  (over  which  the  random 


fluctuations  are  to  be  superimposed)  is  of  utmost  importance. 

At  still  higher  frequencies  (say  X-band  and  beyond)  most  scat- 
terers  appear  to  be  large-scale.  Coherent  background  effects  no 
longer  exist  and  all  scattering  centers  must  be  regarded  as 
belonging  to  the  "fluctuation"  spectrum. 

Figure  1  describes  the  physical  configuration  consistently 
with  previously  presented  arguments.  The  medium  is  modeled  as  a 
slab  of  average  height  h  (which  may  not  coincide  with  the  average 
height  concluded  from  visual  observations);  a  complex  background 
permittivity  eeff  =  e2  which  is  presumed  constant  (more  precise 
models  may  necessitate  more  complicated  plan-stratified  or  possibly 
anisotropic  backgrounds);  and  permittivity  fluctuations  (e  (r,t)  ) 
which  are  random  functions  of  space  and  time  with  supposedly 
known  two-point  space-time  correlation  functions.  *-s  neces¬ 

sarily  complex  accounting  for  absorption  as  well  as  scattering 
losses.  The  ground  is  taken  to  be  a  perfect  reflector  (a 
reasonable  assumption  for  horizontally  polarized  waves).  As 
mentioned  above,  the  random  nature  of  both  vegetation-air  and 
vegetation-ground  interfaces  is  disregarded  at  present.  The 
forthcoming  analysis  presumes  the  dominance  of  geometric -optical 
contributions  to  both  the  incident  as  well  as  the  backscattered 
radiation.  This  is  not  the  case  under  a  wide  variety  of  circum¬ 
stances  such  as  (a)  antenna  within  the  vegetation  (or  within  few 
wavelengths  above  it),  (b)  range  cell  in  the  shadow  region  of  the 
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geometric-optical  refracted  wave  due  to  hilly  terrain,  etc.  It 
can  be  assumed  that  the  coupling  between  the  interfaces  is  small 
(a  reasonable  assumption  in  view  of  the  prevailing  scattering  and 
absorption  losses).  This  allows  one  to  treat  the  forest  as  a 
half-space,  with  ground  effects  accounted  for  as  the  single  and 
double  reflections  illustrated  in  Figs.  8a,  b.  These  ground 
reflections  render  the  determination  of  the  range-cell's  location 
uncertain.  However,  the  introduced  errors  are  expected  to  be 
minimal  in  the  present  context. 

It  is  anticipated  that  the  solution  of  the  random  scattering 
problem  (i.e.,  the  determination  of  the  analytical  relations 
between  the  presumably  known  forest  parameters  and  the  pertinent 
wave  statistics)  via  "distorted  wave  Born  approximation"  would 
provide  acceptable  accuracy  under  the  prevailing  circumstances 

Contributions  to  the  crossed  polarized  wave  do,  in  general, 
stem  from  both  coherent  background  effects  (reflection,  refraction 
and  diffraction)  as  well  as  random  multiple-scattering  events 
(first  order  description  of  the  cross  polarized  field  is  provide 
by  the  second  order  Born  approximation^^).  Since  it  is  a  single¬ 
scatter  mechanism  which  is  believed  to  dominate,  depolarization 
due  to  random  scattering  is  presumably  insignificant.  Coherent 
refraction  and  diffraction  (giving  rise  to  a  lateral-wave)  are 
likely  to  become  major  contributors  to  the  crossed-polarized  wave. 
However,  the  rather  small  discontinuity  at  the  air-forest  inter- 
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face  suggests  that  the  overall  effect  is  small.  This  is  not  to 
say  that  the  original  polarization  dominates  everywhere.  For 
example,  near  the  ground,  or  at  minima  occurring  due  to  destructive 
interference  of  the  horizontally  polarized  wave,  the  vertical 
crossed-polarized  component  may  dominate  locally  despite  its 
smallness.  The  total  amplitude  of  the  incident  field  will  be 
small  and  as  a  result  the  backscatter  from  these  spatial  regions 
is  negligible. 

Section  B  describes  the  physical  configuration  and  formu¬ 
lates  the  scattering  problem  associated  with  the  random  permit¬ 
tivity  fluctuations  superimposed  on  a  lossy  deterministic  back¬ 
ground  slab.  The  background  problem  is  discussed  in  Section  C, 
starting  from  the  exact  field  solutions.  Asymptotic,  geometric- 
optical  forms  are  derived,  tailored  to  the  specific  parameter 
range;  relevant  to  the  prevailing  experimental  conditions. 

Conditions,  restricting  the  validity  of  the  geometric -optical  filters 
are  discussed  and  a  comparison  with  possible  "lateral "-wave  con¬ 
tributions  is  given.  In  Section  D  (and  in  the  Appendix)  we  derive 
expressions  describing  the  scattered  field  consistently  with  the 
distorted  wave  Born  approximation  and  the  observations  made  in 
Sections  A  and  B.  The  mean  received  power,  its  variance  and  one- 
point  distribution  as  well  as  the  field's  temporal  correlation 
associated  with  the  clutter  return  are  calculated.  The  clutter's 
spectral  characteristics  are  investigated  in  Section  E.  Two 
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models,  one  presuming  the  velocity  field  to  constitute  a  multi¬ 
variate  normal  process  and  another  presuming  the  scatterer's 
motion  to  be  quasi -harmonic  over  limited  time  segments  are 
discussed . 

Certain  relevant  aspects  not  presented  in  this  report 
but  which  can  be  readily  described  within  the  general  framework 
of  the  proposed  model  are  mentioned  below: 

(a)  Situations  in  which  the  geometric-optical  field 
no  longer  constitutes  a  major  contribution.  This  is 
the  case  whenever  the  transmitting  and/or  receiving 
antennas  are  situated  within  (or  slightly  above)  the 
vegetation.  A  similar  situation  arises  whenever  the 
range-cell  is  situated  in  the  shadow  region  of  the 
refracted  wave,  owing,  say,  to  hilly  terrain.  The 
"lateral"  wave  constitutes  the  major  channel. 

(b)  Distribution  laws  of  the  scattered  field. 
Amplitude  (power)  fluctuations  and  related  two- 
point  statistics  associated  with  target  and  clutter 
returns . 

(c)  Effects  of  vegetation  (random  slab)  on  target 
spectrum. 

(d)  Smoothing  effects  owing  to  finite  aperture  of  the 
receiving  antenna  and/or  finite  extent  the  target. 
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(e)  Generalization  to  bi-static  situations. 

(f)  Pulse-shape  distortion  and  the  resulting  range¬ 
cell  diffusion. 

The  analytical  details  of  the  model  are,  at  times,  tailored 
to  the  specific  parameter  range  of  interest.  Many  of  these 
restrictions  may  be  waived  if  the  need  arises  under  varying 
circumstances . 

B.  FORMULATION  OF  THE  PROBLEM 

The  pertinent  physical  configuration  is  described  in 
Fig.  1.  Heuristic  arguments  in  support  of  the  random  slab  model 
have  been  presented  in  a  previous  interim  report  and  the 
Introduction  above.  It  should  be  noted  that  the  permittivity 
£2  +  e  (r  ,  t)  is  generally  complex  reflecting  absorption  as  well 
as  scattering  losses  within  the  slab  region.  Since  the 
absorption  and  more  so  the  random  scattering  constitute  frequency 
sensitive  phenomena  it  could  be  anticipated  that  both  an<^ 
e  (r  ,  t)  (and  its  characterizing  statistics)  are  frequency- 
dependent  (dispersive)  as  well.  Also,  the  assumption  of  perfect 
ground  conductivity  while  generally  valid  for  the  horizontal 
polarization  is  unsatisfactory  for  the  vertical  polarization. 

The  propagation  and  scattering  characteristics  of  the  EM 
waves  are  determined  by  Maxwell's  equations: 
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Fig.  1.  Physical  configuration. 
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4i  ,  ^  a|9  a 

X  H-,  -  en  -rr-  +  J-,  v  X  H2  =  e2  ^  (e  E?)  +  J?  (1) 


ll  *1  at  T  ^1 


V  x  E  x  = 


a£. 

"  %  ~ 


V  X  E  2  = 


-  d 


a£2 

o  at 


(2) 


together  with  the  continuity  conditions 


/\  S' 


/X 


-It  — 2t  5  -It  — 2t 


(3) 


at  the  interface  and  a  causality  condition  which  in  the  steady- 
state  case  becomes  the  radiation-condition  at  r-»  =>  . 

It  is  presumed  that  the  temporal  medium  variations, 
characterized  by  the  correlation  time  tc  ,  are  slow  compared  to 
these  of  the  wave,  i.e.,  udqtc  >)  1  where  is  the  angular 
frequency  of  the  presumedly  monochromatic  transmitter. 

Let , 


El  (r,t)  =  Re[EL  (r,t)  e^o*] 
2  2 


(4) 


Hx  (r,t)  =  Rejj^  (r,t)  ejtUotJ 
2  2 


(5) 


(r,t)  =  Re'  J-^  ^  ejuu0t  j 
2  2 


(6) 
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(Note:  since  the  source  is  taken  to  be  monochromatic  J-^  2 

is  time  independent.)  The  substitution  of  Eqs .  (4-6)  into 

Maxwell's  equations  (1,2)  results  in: 


V  x  Ej. 


V  x  Hx 


V  X  ^2 


V  X  H2 


-1  uu  n 

J  00 


(l  +  Ji;  If)  Hi  +  Ii  w 
1  Vl  (L  +  TE)  h 

(“it  *)■’ 

(*  +  If)  te2  h  +  e  £2]  +  —2  <r> 


00 


(7) 

(8) 

(9) 

(10) 


from  which  the  following  vector  wave  equations  are  readily  derived: 


V  x  V  x  ex  -kj^l  +  at)  -1  ju)oMo  ^1 

(2 

1  +  1%  ^2  =  j%^o  —2  +  %Uo 


(11) 


(12) 


2  2 

where  k.  =  a)  u  e- 

1  o  o  i 

2  2 

Terms  associated  with  the  slow  temporal  derivative  may 
be  safely  omitted  due  to  the  smallness  of  (vc) "  “  I®  << 1  • 

(here  Auu  denotes  the  frequency  deviation  of  the  wave  from  UJQ)  • 
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2  d 

For  example ,  comparing  the  terms  | e  E2  |  with  'jir  st  (e  —2^ '  one 
has 


<e  1  'e  ^  (ff  ^2  +  6  I2' 


2  /Be 


BE, 


(w  T  )‘1  +  M  <(  1 

V  O  C7  UD 

O 

consequently,  we  may  take 


(13) 


1  +  «  1 
ju)o  at 


and  equations  (11,  12)  reduce  to 


V  x  v  x  E-,  -  k-.  En  =  - j  uj  u  J-, 

— i  i  — i  00  — 1 


(14) 


V  x  V  x  E2  -  E2  =  -ju>0U0  J2  +  uu^po  e  E2  ,  e  =  e  (r,t)  (15) 

in  which  time  appears  as  a  parameter. 

Let  2  denote  the  Dyadic  Green's  function  defined  by 
the  equations, 


V  x  V  x  Gx  =  0 


v  x  v  x  G2  ■  k2  G2  =  15  (r - r ' )  -h(z(0 

the  radiation  condition  at  infinite  and  the  appropriate 
continuity  conditions  at  z  =  0.  Here,  I  denotes  the  unit  dyadic, 


z>0 

(16) 

-h(z(0 

(17) 

12 


denote 


Let  E^2  =  2  (®=0)  and 


=  E  -E(°) 
2s  — 1 , 2  — 1 , 2 


the  unperturbed  (background)  and  the  scattered  waves,  respectiv 
From  the  basic  definitions  and  Eqs.  (14,  15)  one  obtains: 


’  *  ’  *  hs  -  ki  h*  *  0 


v  X  vxE2jS  E2s  -  po  e  (r  ,t)  E2 

which  may  via  Eqs.  (16,  17)  be  converted  to  the  convenient 
integral  forms 

Eis  "  *20  ,oJ(d3  rl  *<El’t>  Si  <£.  Ei>  ’  —2  (Ei’t) 

—2s  =  %  doXsd3  rl  e(-l5t)  S2  <E»  Ei>  '  E2  (Ei»t) 

Equations  (20,21)  constitute  integral  equations  which  generally 
cannot  be  solved.  The  simplest  applicable  approximation  to 
these  equations  is  the  so  called  "distorted  wave  Born  approxima¬ 
tion"  in  accordance  with  which  one  replaces  EL,  (r^,t)  in  the 
integrands  of  Eqs.  (20,21)  with  E^^  (r.pt).  The  quantities 

and  G-^  2  are  associated  with  the  background  problem  which 
is  discussed  in  the  next  section.  It  should  be  noted  that  the 
validity  range  of  the  Born  approximation  (in  absence  of  any 
wave  distorting  background)  has  been  thoroughly  iiv  estigated 
(e.g.,  ref.  (2))  and  further  elaboration  seems  redundant. 


ely . 

(18) 

(19) 


(20) 


(21) 
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C.  THE  BACKGROUND  PROBLEM  (EVALUATION  OF  G  AND  E^) 


For  the  sake  of  simplicity  the  ground  effects  are 
presently  ignored;  they  will,  however,  be  added  separately  at 
the  end  of  this  section.  The  problem  to  be  considered  is  that 
of  a  dielectric  half-space  excited  by  an  arbitrary  dipole.  It 
is  assumed  for  definiteness  that  the  source  is  located  in  region 
2  (the  half-space  z)0).  To  obtain  solutions  applicable  for  a 
source  situated  in  region  1,  one  simply  replaces  1  by  2,  2  by  1 
and  z  by  -z .  The  excitation  properties  of  a  dielectric  half¬ 
space  by  a  dipole  are  well  understood,  extensively  elaborated, 
upon  in  the  literature  and  are  regarded  as  textbook  material. 
Despite  their  simplicity,  no  exact  closed-form  solutions  are 
available  and  one  must  resort  to  an  asymptotic  description  which 
is  readily  interpretable  in  geometric-optical  terms. 


1. 


The  Exact  Solution. 


The  satisfaction  of  the  equations 


(22) 


with  the  continuity  conditions 


(23) 


and  the  radiation  condition  at  infinity  is  required. 
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The  boundary  value  problem  is  conveniently  describable  by 
the  vector  potential 


H-^  =  v  x  A-^  ,  E-^  = 
2  2  2 


■JVo 


1  + 


V  V 


^1 

2 


-1 


with  A^  obeying  the  reduced  wave  equation: 
2 


(v2  +  k2)  Ax  =  0 

(v2  +  k2)  A2  =  -J2  6 (r  -  r ' ) 


Assume  the  following  plane-wave  representation. 


Al(r)  =  fa  kt  h  <kt’z)  e’j“t  *  ~  ~  ) 

1  C2TT)2  J  J  c 


hence 


[fr  +  “i]  -i  ’ 0 

-2  "  -2  6<Z  ‘  2,)’  *}  ‘  kJ 


(24) 


(25) 

(26) 


(27) 


(28) 


(29) 


where,  £  =  i:  -  z^  z  >  2.'  ~  ~^oz '  anc^  — t  denotes  the  transverse 
wavevector . 


A^  are  of  the  form 
2 
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£2  =  j"x  e-j*2fe-z')+  r  ej*2  <Z+Z'>J. 


^0 


(30) 


Ai  =  T  •  A 
—1  ~  — o 


,-j* iz  +  j*2z ' 


(31) 


where 


r  =  th 


+  Te 


z  z  -  — 7T  z  k 
— o  — o  ,2  — o  — t 

Kt 


T  =  T 


H 


X1 

i„  +  -4 

z  kj 

+  T„ 

z 

z 

x2  1 
•  4  z  k. 

-t  k2 

-CJ 

1  E 

— o 

— o 

4  _tJ 

(32) 


(33) 


consistently  with  the  continuity  conditions  at  z  =0,  and 

A  =p- 

-O  2jx2 

consistently  with  the  prescribed  dipole  excitation.  I  represents 
the  transverse  dyadic 


(34) 


I,.  =  I  -  z  z 
^  — o  — o 


(35) 


Fr  (Th)  and  r£  (Tg)  denote  the  reflection  (transmission)  coefficients 
associated  with  the  H  and  E-mode  constituents,  respectively. 

They  are  given  by: 

2x„ 


x2  *  X1 
H  x0  +  x. 


TH  Xo  +  X- 


(36) 
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rE 


x 

kl  X2 


i?"x 

K1  X2 


T„  = 


2k*  x2 

~1 - 1 - 

kl  x2  +  k2  X1 


(37) 


It  is  customary,  although  for  some  purposes,  inconvenient, 
to  make  a  transition  from  the  plane-wave  representation  (27) 
which  involves  a  two-fold  integration  to  a  cylindrical  represen¬ 
tation  involving  a  single  integral.  Omitting  intermediate 
algebraic  steps,  one  obtains: 


z-z 


—2 (— ’  £'>  -  SSJ  /ds  ^  "02>  (P|£'£,|>  [l  e'J*2 
-00 

+  r  (p)ejx2(^/)]  .  4  +|a  vty*dp^H<: 
X  0|p-p'|)  r2  (P) 


(2) 


-00 


.  J2 


(38) 


Ai(r,r')  -  jij  J  dB  (B|z-z'|)  Xl  (B)  .  j2 


Z 

— o 
3F 


y*dp  ^  ho(2)  t2  (p) 


x  e 


-j  xiz+j  X2Z  / 


^2 


(39) 


where  =  7  •  ^  denotes  the  transverse  gradient,  P=|kt|  and 
~1’  ~1’  ^2  and  ^2  are  defined  by  the  relations: 
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(40) 


T.  =  T„  L  +  Tr  z  z 
~1  H  ~t  E  — o  — o 


~L  “  rH  ~t  +  rE 

T2  "  7  (*1  TH  -  H2  TE> 

r2  =  "  (rH  +  rE) 


(41) 

(42) 

(43) 


The  above  quantities  are  related  to  the  previously  defined 
F  and  T  via 

E  -  Jj.  +  r2  ^  kt 


T  =  T-,  +  T0  z  k. 

~  ~1  /.  — o  — t 

Equations  (38,  39)  cannot  be  further  reduced  and  as 

expected  the  resort  to  approximation  is  unavoidable. 

2.  The  Asymptotic  Evaluation.  The  Geometric -Optical 

Contribution. 

We  first  replace  H  ^  (3  ljL~.fi/ 1)  in  Eqs .  (38,  39)  by  its 
asymptotic  values 

ho<2>(x)^  e~jx+j^  [1+0  (£)  ] 

resulting  in 

A2  (r,  r')~ 


(44) 

(45) 


eJ5  JS 


8rrj 


vTsFl  I 


/ 


de  r1 

X2 


I  e‘jPl  +  r,  e‘jP2 


+  z  p 
— o 


+» 

/ 


dp 


1 


3/2 


r2  (P)  e 


-jP2 


—2 


(46) 
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(47) 


where, 


P]_  (P)  =  P|£-£;  |  +  *2  \z~z'\  (48) 

P2  (p)  =  P  I  1  "  H2(z-te/)  (49) 

P3  (P)  =  P|£-£7  |  +  z-x^z'  (50) 

and  the  asymptotic  relation  -j  pp^  =  |  <  |  is  a  unit 

vector  in  the  £.-£_*  direction)  has  been  utilized. 

The  major  asymptotic  contributions  to  the  p-integrals  (k  the 
large  parameter)  come  either  from  the  vicinity  of  the  respective 
saddle-points  (geometric  optical  contributions)  or  (within  some 
restrictive  spatial  regions)  from  the  appropriate  branch-points 
(lateral  waves).  Situations  under  which  lateral-wave  contributions 
are  expected  to  dominate  the  channel  (see  Introduction)  are  omitted 
and  will  be  reported  separately.  Their  relative  significance, 
however,  is  discussed  at  the  end  of  this  Section, 
a)  The  Saddle-Point  Contributions 

Consider  the  integral 

-jP(P) 


I  = 


/ 


dp  F(p)  e 


(51) 
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assuming  the  existence  of  a  single  first  order  saddle-point  one 
obtains 


cp 


V|P"^S> 


-jP(Ps)  -j l  sgn  (P"(3g)) 


f(Ps)  e 


(52) 


where  the  saddle-point,  3  is  determined  by  the  saddle-point 

o 

condition 

p'CSg)  =  0 


(53) 


and  the  prime  denotes  the  derivative  with  respect  to  the  argument. 
One  can  readily  apply  the  asymptotic  form  (52)  to  Eqs .  (46-47). 

The  results  are  easily  interpretable  in  geometric -optical  terms 
and  could  alternatively  be  derived  via  a  direct  use  of  geometric- 
optical  arguments.  These  are, 


-j 

4(1)  - 


-j  1<2  (L^  +  Lr) 

=  VdY  -  Lr)  [El  <V  +  2o  £o  r2  <».)J  •  h 


P 


(t)  _  e 


-j(k2Li  +  k1Lt)  rgg- 


4tt  Li 


V1F  [ll  +2o£o  es  T2(0s)J 


(54) 

(55) 


—2 


(56) 


The  length  segments  L^,  Lr,  and  Lt  are  shown  in  Figures  (2,  3); 

da  and  da  denote  differential  cross  sections  of  the  transmitted 
o 

ray  tube  at  the  observation  point  and  at  the  interface,  respectively: 

af  =  (l+T1  r1  Wl+r*  (57) 
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Fig.  2.  The  geometric  -  optical  reflected  wave. 


Fig.  3.  The  geometric  -  optical  transmitted  wave. 


21 


where  the  relationship 

{3  =  k0  sin  9 .  =  k,  sin  0 

S  Z  1  1  t 

has  been  used,  and  from  Eqs .  (36,37) 


(58) 


cos  9^  -  \^2  cos  ^ 

*E  ^s^  "  cos  0^  +  k2  cos  0^*  TE  ^s^  “  k^  cos  et  +  k2  cos  9(. 

(59) 


2  k,  cos  0. 


k0  cos  0.  -  k,  cos  9, 


^2  w  tv  ^  vuo  ^  £. 

(gs)  -  k2  cos  0^  +  k^  cos  0^"’ 


2  k2  cos  0^ 


^H^s ^  k^- cos  0^  +  k-^  cos  0^ 

(60) 

The  geometric -optical  electric-field  is  derived  via  Eq.  (24). 

Neglecting  rapidly  decaying  field  constituents  one  obtains  the 
incident  electric  field. 


■(i) 


-j 


o  "o  4tt 


(I  -  r  r  ) 
~  —o—o 


—2 


(61) 


where 


r  -  r 


-j,pl  <6s>  *  -Jk2  £o>  £o  *  Tr  ~-~r ' 

the  reflected  electric -field, 

-jk2  (Li  +  Lr) 


E 


-JU)_  p, 


o  ^o  4tt  (L^  +  Lr) 


a  -  ££,) 


r  <pH) 


—2 


(62) 


(63) 


where 


jvp. 


Os)  =  -jk 


/\ 

2  ^o’ 


/\ 

r 

— o 


£-£*  +  (z+z  ' ) 
l£~£J  +  Zq  (z+z  ')  | 


(64) 
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/\  H2  ^e) 

r  (S.)  ■=  rH  (Ss)  it  +  t-e  (Ss)  ^  -  A-S- 


rH  <es>  +  rE  <Ps>  5o£o 


(65) 


and  finally  the  transmitted  electric  field, 

-j  (k2Li+k1Lt)  r 
- jW  VI  - - j—r: -  \ 

J  o  o  4tt  L_.  i 


rar 


-sp 


“3a  C1  ’  I0Eo) 


I(Pe) 


—2 


(66) 


where 


and 


7~  "j vp3  <Bs>  “  -Jkl  V  Eo  ‘  £o  sln  ®c  +  So  cos  ec 

i  <*>s>  *  th  <f>s>  it  +  TE  (Ss>  2oZo 

+  ^  [*1  <Ss>  TH  ^s>  -  *2  <Ss)  TE  (Bs>]  5o£o 


(67) 


(68) 


r^,  r^  and  r^  geometrically  represent  unit  vectors  along  the 
incident,  reflected  and  refracted  rays.  The  dyadics  1  -  r^r^, 

1  -  r^r^  and  1  -  F^r^,  represent  transverse -project ion  operators 
with  respect  to  the  corresponding  rays. 

With  the  geometric -optical  properties  of  the  refracted 
wave  (66)  available,  the  desired  information  concerning  both  the 
incident  wave  and  the  dyadic  Green's  function  is  readily  deduced. 
For  reasons  of  convenience,  notational  changes  are  introduced  as 
indicated  in  Figs.  (4,  5).  The  geometric-optical  Green's  function 
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Fig.  4.  Incident  wave. 


Fig.  5.  Green’s  function. 
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and  the  incident  wave  (excited  by  an  arbitrary  dipole)  are  givqn  by: 
G 


”j  (k^L^+k^l^ )  f 

h  (£-£i>  ~  5 — &nr2 - V 


da  <2> 
o 

m 


da 


wi> 


-12 


(69) 


-3  (k1L1+k2L2)  /da(l) 

:  (£x>  ~  "3  Vo  TZTl{ - \  ^ITJ  (it  " 


where 


Sl2  " 


~21  ‘  -1 


(12)  _  (12) 

t<12>  I  +  Tp12^  zz  +  — — S-  2  E 

H  — t  E  -o-o 


(70) 


T 


Vo 


e=e 


s 

(71) 


-21 


th21)  It  +  te21)  z<A>  - 


*  T(21>  -  H  T<21> 

2  iH  H1  LE 


Vo 


0=0 


and 


(12,21)  _  2w2,1  t(12,21)  .  2kl,2  *2,1 

H  "  x,  +  x  '  E  '  k2  x„  +  k2  x, 


s 

(72) 


(73) 


ll T  2 


"1  ^2  2  1 


The  radius  vectors  r  and  r^  denote  the  location  of  the  source 
(generating  the  incident  wave)  and  the  scattering  point,  respectively. 
The  ray-tube  cross  section  ratios  are  now  given  by: 

da^-2) 


da 


_  1 1  ,  L2 , 1  kl,2 

ttjj  ~  l1  + 


i  + 


L2 , 1  kl,2  COS^l,2 


J1 , 2  k2 , 1  cos^e. 


(74) 


-0  v  -  2,1 

Further  simplifying  features  which  stem  from  the  presumed 
smallness  of  the  angle  6  =  ^  -  6^,  are  discussed  below. 
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Specifically,  the  validity  of  the  inequality 


sin  6«o  «  1  - 


rf 2  «  1 


(75) 


is  assumed.  Here,  p  =  k2/k^  denotes  the  index  of  refraction  of 
the  half-space  z  <  0  .  (Typically,  ri  «  1.05  and  (75)  becomes 
6  «  0.3  radians  «  17—,  often  a  reasonable  restriction  for  ground 
radars).  Under  the  restriction  (75)  one  may  readily  verify  that 


2 

cos  90  i  -2 _ 2  .  ,  - 

2  _  1  -  p  cos  o  ..  1  -  T| 

t: - rr_,2“ - 2 


»  1 


cos  9, 


sm 


sin  6 


(76) 


consequently,  the  area-mappings  (74)  simplify  to: 


da(1)  _  I  +  L2  kl  „  i 

—nr  ~ 1  +  it  E7  "  1 


da 

da 


o 

(2) 


da 


2 

L,  k~  cos  9« 

si  ^7 


L1  ^2\  /  ^ i  ^2  cos 

L j,  k^  J  I  k^  cos  9-j 


(77) 


(78) 


where  the  last  of  the  approximate  forms  in  Eqs.  (77,78)  presumed 
the  largeness  of  L-^  with  respect  to  L2 .  More  caution  must  be 
exercised  in  approximating  the  phase  function: 


P  =  k1  Lx  +  k2  L2  (79) 

The  following  geometrical  identities  hold: 

|  £~£^  I  =  L-^  sin  9^  +  L2  sin  92  ,  £  =  r  -  zz^,  p^r^-z-^z^  (80) 
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z  =  L-^  cos  9^ 


(81) 


"Z1  =  L2  COS  e2 

d02  k.  cos  9. 

From  Eq.  (76)  it  follows  that  =  k2  cos  ~2  <<:  l-  hence'  £or 

large  |£-£jJ  (small  6)  the  angle  02  stays  virtually  constant 
approaching  the  critical  angle  of  refraction: 

^1  1 
sin  e2  Miq  -  $ 

It  follows  that : 

-z. 


(82) 


(83) 


"i  ”2 


(84) 


To  determine  we  substitute  Eq.  (81)  into  (80)  resulting  in 

k  k  1/2 

l2-£ll  “  +  l2  sin  91  -  <L1  +  h  [1  '  1  <85> 


upon  expansion  of  the  square  root  one  obtains 

,2 


klLl  “  kl  I  £-“£-1 1 


i  +  (tr) 


vH+  ... 


k2 
1  T 

q  L2 


(86) 


Neglegence  of  the  fourth  power  contribution  to  the  phase  function 
is  justified  subject  to  the  constraint 


(Li) 


«  1 


which  is  an  acceptable  assumption  under  the  prevailing 


(87) 


-3, 


circumstances.  (Typically,  at  UHF:  k^!£-£^|(^— )  <  0  (10  ) 
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From  Eqs.  (85,87)  it  follows  that 

l£-£ll  -  a,  +  i  V  [l  -  \ 


2 


=  L-, 


1 

2  L, 


+  ^l2  -  7S  Vq;>' 


and  presuming 


(88) 


1  ^2  L2  <<: 

(once  again  a  reasonable  assumption  under  typical  conditions  at 
UHF: 


(89) 


1 

1 


(f") 

L1 


^  0 


do"2) 


one  readily  solves  the  resulting 


(90) 


The  square-root  expansion  in  (90)  brings  about  a  constraint 
already  imposed  by  (87). 
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We  finally  have 


L1  =  "if  +  1 


TT  *  lfi.-fi.il  "  “ 


l£-£l 


+  i  z‘ 


2  l£-£x 


(91) 


z2L, 


1  2 

The  last  form  presumes  the  smallness  of  the  term  -s-ki  - * 

2  1  i£-£li2 

a  restriction  equivalent  to  (89).  The  phase  function  is  given  by: 


P  ~  k 


1  +  ?  1#^t] 


r^7  z. 


The  last  result  is  simply  interpretable  as  illustrated  in  Fig 
Equation  (79)  may  be  rewritten  as: 


P  =  kx  yj  +  z2  +  L x  "  V  !£-£!  | 2  +  z2 


4"  1^2  1 2 


The  approximation  (92)  amounts  to  replacing  the  bracketed  term 
by  -L.2  sin  ©2  (see  Fig.  6)  and  to  the  subsequent  expansion  of 
the  square  root,  retaining  the  two  leading  terms.  Equation  (92) 
follows  directly. 


29 


(92) 


6. 


(93) 


Fig.  6.  Interpretation  of  the  phase  function. 
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The  polarization  properties  of  the  incident  wave  and  the 
Green's  functions  are  determined  by  the  operators  (1  -  r^r^)  *  %21 

and  (1  -  r^r^)  •  T^2  (see  Eqs.  (69-72)),  respectively.  These 
operators  simplify  substantially  in  view  of  (75).  One  has, 

(i-lollol)  ‘  T12  =[l  -  sin2  0j_  -5o5o  cos2  ei  (£0^0  +  SoLo) 


sin  0 


1  cos  9i]  •  I12  “  [i  -  £<A)]  • 


T12  “  TlT'  (it  Vo) 


(12) 


+  T 


(12) 


[ 


z  z 
-0-0 


cot  0 


(94) 


where  terms  of  order  cos  0^  =  sin  6»6  have  been  omitted.  The 
depolarization  may  also  be  neglected  owing  to  the  smallness  of 


cot  $2  ^  v  r|  -1  «  1 
Similarly, 

(A  -  *  ^21  =  [i"AA>  sin2e2'^o  c°s202-(£ozo  +  z^) 

sin  02  cos  02]  •  T,,  w  t^21)  _  £^) 

+  T(21)  z  z 

+  le 

2 

where  once  again  the  smallness  of  the  parameters  6  and  1  —  rj  was 
presumed.  With  restriction  (76)  borne  in  mind,  the  following 
simplified  forms  of  the  transmission  coefficients  results: 

T^12^  »  T^12^  *  2 


(95) 


(96) 


31 


(97) 


T(21)  *  1  t(21>  «  2  —  6l  =  2  -z. 

H  7  E  11  cos  9 2  V7-T  L1 

The  substitution  of  Eqs.  (77,78,92,96,97)  into  Eqs.  (69,70), 
results  in  the  following  simplified  geometric -optical  expressions: 


Si  (2L>2Li) 


z[l  ~ 

2  W  n2  -  i  L2 


2 

z 


l£-£ll 


(98) 


and 


(— 1>  j^o^o 


zCl  -  £qPq1  *  ±l 

2tt72-11 


e-Jki  [l£-£il 


+ 


1 

7 


1  £"£i  I 


(99) 


where  a  decay  term  exp  (-otl^)  has  been  written  explicitly  (not 
as  the  imaginary  part  of  k^) ,  and  ^  are  given  by  Eqs.  (84,91). 
(The  last  two  terms  in  Eq .  (91)  may  be  totally  ignored  in  the 

amplitude  factor  of  Eq.  (99).) 

For  a  horizontal  dipole,  we  have 


(1  -  £o£o)  ‘  li  =  V^o  *  ±1  =  SBq  Jlcp  =  J1  cos  ^ 


(100) 


(see  Fig.  7) 

The  incident  field  may  be  further  generalized  if  one  cares 
to  account  for  the  finite  aperture  of  the  antenna.  Equation  (99) 
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=  (AS,)h 


CENTER 

BEAM 
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may  be  rewritten  as : 


,(o) 


2z 


■Jki 


'£"£l'  +  ?  Ti^T  '^r~lzi 


-aL„ 


2o 

(101) 


where  F(cp)  is  simply  related  to  the  incident  field  as  follows: 
F  (cp)  =  F(9,cp)0=o 

E(i)  ~  F (9 ,cp) 

From  basic  definition  of  the  antenna  gain  (ignoring  medium 
effects)  via, 


(102) 

(103) 


g(e 


,») .  ieg* 


sL 


(104) 


t  o  t 

It  follows  from  Eqs.  (102,  103,  104)  that 

I f (cp)  | 2  «  Pt[g(e,cp)]e=0 


(105) 


where  denotes  the  total  power  emitted  from  the  antenna.  It 
should  be  recalled  that  the  antenna's  finite  aperture  not  only 
plays  a  role  by  controlling  the  incident  radiation,  but  also  by 
receiving  the  backscattered  signal.  The  necessary  related 
modifications  are  described  in  the  next  section. 
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b. 


Ground  Effects 


The  presence  of  the  ground,  so  far  totally  ignored. 


can  be  accounted  for  without  difficulty.  For  the  horizontally 
polarized  component  the  ground  could,  quite  accurately,  be 
modeled  as  a  perfectly  reflecting  surface.  Furthermore,  multiple 
reflections  between  the  two  interfaces  may,  safely,  be  neglected 
owing  to  the  lossy  character  of  the  vegetation  slab.  A  single 
ground  reflection  of  either  the  incident  wave  or  the  Green's 
function  (or  both)  is  expected  to  constitute  an  adequate 
description  of  the  phenomenon.  Geometric -optical  contributions 
to  be  dealt  with  are  illustrated  in  Figs.  (8a  and  8b).  The 
ground  reflected  constituents  may  be  written  by  inspection  from 
Eqs.  (98,99,  or  101)  by  a  proper  phase  inversion  (for 
horizontal  polarization)  and  by  replacing  -z-^  with  2h  +  z^. 

We  have. 


(106) 


(Note:  the  z^^  element  does  not  experience  a  phase  reversal 

but  it  plays  no  significant  role  in  future  calculations.) 
Similarly, 
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(a) 


(b) 


Fig.  8.  Geometric  -  optical  contributions. 
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Er0>  (£!>--£„  F(op) 


2z 


>-Jkl 


V  n2-l  L? 


JlP-P! 


I  + 


l£"£ll 


1 


—* —  "1  a  (2h  + 

rf-1  (2h  +  z,)  - 

J  n Z^T 


zl> 


(107) 


c.  Comparison  with  the  Lateral  Wave  Contribution 

Under  certain  conditions  the  geometric-optical  field 
(Eqs .  (98,99))  does  not  represent  an  adequate  description  of  the 
incident  radiation  or  the  Green  function.  Specifically,  as  the 
radiating  source  or  the  point  of  observation  approach  the 
interface,  i.e.,  as  z -»  0,  the  geometric-optical  field  vanishes 
and  clearly,  it  no  longer  constitutes  the  major  filed  contribution. 
Whenever  the  transmitting  and/or  receiving  antennas  are  situatec 
in  the  proximity  of  or  below  the  (air -vegetation)  interface,  it 
is  a  diffracted  contribution,  the  lateral  wave,  which  dominates. 
While  work  has  been  done  to  account  for  such  situations,  its  report 
is  deferred  for  the  future.  Under  the  conditions  stated  in  (75) 
the  dipole -excited  lateral  wave  is  given  approximately  by: 
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where  LT  denotes  the  lateral  path.  Consequently,  comparison  with  eq.  (99) 
L 


yields, 


1  1 
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X 


(109") 


Typically,  for  antennas  situated  several  wavelengths  above  the  vegetation 

air  interface,  the  lateral-wave  may  be  safely  neglected. 

D.  THE  DISTORTED  WAVE  BORN  APPROXIMATION.  THE  RECEIVED 
CLUTTER  POWER  (MEAN  AND  VARIANCE) 

The  title  ’’Distorted  Wave  Born”  generally  refers  to  a  class  of  single¬ 
scattering  approximations  in  which  the  propagation  features  of  both  the  in¬ 
cident  wave  (on  the  way  to  the  scattering  volume)  as  well  as  the  scattered 
wave  (on  its  way  to  the  receiver)  are  modified  consistently  with  a  suitably 
defined  background  (in  our  case  a  lossy,  uniform  dielectric  slab).  In  the 
following  analysis  we  refrain  from  generalizing  any  of  the  results  beyond 
restrictions  applicable  for  the  specific  configuration.  Several  such  generali¬ 
zations  however  are  within  easy  reach.  From  eq.  (21)  together  with  (98,  101) 

one  readily  obtains: 

Similar  observations  were  made  by  T.  Tamir  in  an  ECOM  Memorandum 
(Summer  1970)  entitled,  ’’The  Electromagnetic  Field  Radiated  Above  a 
Forest  by  an  Antenna  Embedded  in  Vegetation". 
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where  L^  has  been  replaced  by  in  the  amplitude  term  (see  eq.  91  and 
Fig.  7). 

Let  us  focus  our  attention  on  the  scattering  contributions  from  a 
volume  i  V.  (Fig.  7)  such  that  PQAcp^  is  large  compared  to  the  transverse 
correlation  distance  (but  small  enough  to  obey  restrictions  (113,  114).  We 
have, 
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together  with  the  assumptions 
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result  in, 


A£si(r,  t) 


a)  u  z 
o  *0 

T 2  “  4 

TT  (  T|  -  1)  Pq 


e-%  Po  •j2klPo  F«Pi>  2oi 


f d3rj  e(£i+r1,t)ej2kl-&i  '  Jl  +  j2*l  z,  + 


2  az. 


AV- 


'/r~rr 


(1151 


where  it  has  been  furthermore  assumed  that  cgQ  and  F(cp)  stay  virtually  con¬ 
stant  over  the  volume  of  integration  (A  V.). 

Before  proceeding  with  the  analysis  let  us  generalize  eqs.  (110  and  115) 
to  include  the  ground  effects.  We  have  (via  eq.  20), 
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where  G  j,  G ^ R,  and  are  given  in  eqs.  (98,  101,  106  and  107). 


Explicitly: 
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The  effects  of  the  perfectly  conducting  ground  are  totally  contained  in  the 
bracketed  term  (to  be  denoted  y(z^))  its  replacement  by  unity  leads  back  to 
eq.  (110).  Analogously  to  eq.  ( 1 1  5)  one  now  has , 
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It  is  important  to  point  out  an  error  introduced  in  writing  eq.  ( 1 1 6 ) .  The 
scattering  volume  Vg  is  determined  (longitudinally)  by  the  radar's  gating. 

In  eq.  ( 1 1 6 )  we  identify  four  distinct  and  readily  interpretable  terms,  each 
associated  with  a  somewhat  different  range  cell.  The  forms  (116-118)  totally 
ignore  such  differences. 

Since  the  volume  AV.  has  been  selected  consistently  with  the  requirement 
that  PQAcp^  be  large  with  respect  to  the  correlation  length,  it  follows  th^t  the 
scattering  volume  AW  and  AV.  (j^i)  are  uncorrelated,  hence 


<  A  E  .  -  A  E  .  >=  <  Ia  E  .1  >6. 
— s  l  —  s  j  ^  — s  1  1  1 


(119) 


where  <  >  denotes  ensemble  averaging.  We  now  compute  the  field’s  temporal 
correlation  from  eq.  (118) 
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where, 


c(£l*  Lz’  T)  =  <  e(£1*  6  (£.2’  fc  +  T  )  > 


(120) 


(121) 


is  the  space-time  correlation  of  the  presumedly  stationary  random  process  e. 
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It  is  convenient  to  separate  the  vertical  from  the  horizontal  coordinates 
in  eq.  (120)  and  transform  to  the  center  of  mass  system.  This  results  in 
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where  ^  =  p  ^  -  p^,  p  =  (p^  +  p^)  and  AS^  =  PQdACfK.  In  the  transition  from 

eq.  (120)  to  (122)  it  has  been  assumed  that  the  linear  dimensions  of  AS^  (i.  e.  , 
d,  the  range  cell  and  PQAcp^  are  large  compared  to  the  transverse  correlation 
length.  It  has  been  further  assumed  that  the  process  e  (_r,  t)  is  statistically 

homogeneous  and  isotropic  in  the  transverse  direction.  The  integration  over 

2  2a 

d  p  is  evaluated  trivially  resulting  in  AS^,  while  integration  d  p  (for  large  AS^) 

is  recognized  as  the  Fourier  transform  of  C  with  respect  to  £  (i.  e.  the  trans¬ 
verse  spectral  density  function)  to  be  denoted  by 
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Owing  to  the  presumed  transverse  isotropy  r  becomes  an  exclusive  function  of 
h=  \k\9  and  as  is  clear  from  eq.  (122)  it  is  the  value  of  $  at  h  =  |2k^  |  =  2k^ , 

which  is  of  interest. 

The  mean  received  power  (<P^>)  can  now  be  calculated.  From  eq.  (105) 
and  the  presumption  of  a  matched  antenna  one  has, 
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where  <S.>  denotes  the  mean  Poynting  vector  at  the  receiving  antenna;  inte¬ 


gration  over  the  pattern  yields 
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V  1  -  n 

Equation  (126)  may  be  further  simplified  if  one  is  willing  to  except  the  as  sump- 

/S  /\ 

tions  that  the  medium  is  statistically  homogeneous  in  z  (i.  e.  ,  $  =  $  (2k^,z^-z2,  t)), 
that  the  corresponding  correlation  length  is  small  compared  to  the  forest’s 


effective  height  (h)  and  that  ground  effects  are  negligible.  Arguments  analogous 
to  these  used  for  the  transverse  domain  lead  to: 
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where, 
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It  should  be  pointed  out  that  the  assumptions  allowing  the  transition  from  eq. 
(126)  to  eq.  (127)  may  be  seriously  challenged.  However,  the  ensuing  sim¬ 
plifications  are  substantial  indeed. 

Equation  (127)  indicates  that  the  mean  received  power  decreases  as  the 

seventh  power  of  the  range.  Experimental  verification  of  this  dependence 

is  demonstrated  in  Fig.  9  which  is  a  plot  of  received  clutter  power  as  a 

function  of  range  at  UHF  for  a  uniform  pine  forest. 

The  scattered  field  given  in  eq.  (117)  or  (118)  constitutes  (via  the  central 

limit  theorem)  a  normal  process  provided  that  V  (or  AV.),  the  scattering 

volume,  is  dimensionally  large  in  comparison  with  the  correlation  distance 

characterizing  e(r,  t).  It  can  be  shown  (see  Appendix)  that  the  processes 

E.  =  R  (AE  .)  and  E~  =  Im(AE  .)  are  normally  distributed  with  equal  variance 
I  e  si  (L,  — s  i 

Z 

(to  be  denoted  by  a^,)  and  are  uncorrelated.  From  which  it  follows  that  the 
field’s  amplitude  (A.,)  and  phase  are  describable  by  Rayleigh  and  uniform 
distributions,  respectively.  Consequently  (see  Appendix), 
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Fig.  9.  Clutter  power  vs  range. 
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where  <P^>,  the  mean  power,  is  given  in  equations  (126)  or  (127). 

All  information  concerning  the  field's  temporal  correlation  or  its  spectrum 
is  contained  in  the  normalized  function. 


(128) 


where  M  is  a  normalization  constant,  conveniently  selected  in  accord  with  the 
condition  C^o)  =  1. 

hi 

E.  THE  CLUTTER  SPECTRUM 


This  section  is  reserved  to  a  further  investigation  of  C^t)  as  given  by 


eq.  (128)  and  the  associated  spectral  density  under  various  circumstances  to 
be  specified  below. 

The  temporal  fluctuations  associated  with  the  random  process  e(r,t) 
can  be  thought  of  as  caused  by  the  motion  of  distributed  scatterers  characterized 
by  a  random  velocity  field  V  (r,  t).  Thus,  assuming  a  conservative  ensemble 
of  scatterers  (that  is,  assume  that  no  scattering  centers  disappear  from  the 
collection)  one  concludes  that  the  relationship 


(129) 


t 


is  satisfied.  It  is  further  assumed  that  the  processes  e  (r)  and  V  (r,t)  are 
statistically  independent.  This  assumption  is  not  totally  baseless  since  the 


fluctuations  of  e  are  primarily  a  consequence  of  the  spatial  distribution  of  the 


vegetation  (and  therefore  virtually  independent  of  wind  conditions)  while  the 
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nature  of  depends  greatly  upon  prevailing  winds.  As  a  result  the  averaging 
over  the  processes  e  and  V  can  be  carried  sequentially.  One  has, 

C(R,  r,  t)  =  <  e(r,t)  e  (r+R,  t+T  )>  = 
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where  the  process  €  was  presumed  statistically  homogeneous  and  stationary. 

Let  $  denote  the  (spatial)  spectral  density  corresponding  to  C  .  It  follows 
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from  eq.  (130)  that 
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where  V  denotes  the  velocity  component  alongjt  and  the  integration  range 
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t^  s^  t+T  was  shifted  to  o^  s^  t  in  view  of  the  presumed  statistical  stationarity 
of  the  process  V(r,  t).  If  the  process  Y_  is  also  transversely  homogeneous,  the 
term  <  >y  in  eq.  ( 1  3 1 )  becomes  independent  of  the  transversely  homogeneous 
ordinates  (but  is  generally  a  function  of  the  vertical  coordinate).  Separating 
the  transverse  integration  from  the  longitudinal  eq.  ( 1  3 1 )  may  be  written  in 
the  form: 
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Using  the  notation  of  eq.  (123), 
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Owing  to  the  near  grazing  incidence  of  the  wave  we  know  that  the  major  con 
tributions  will  come  from  ranges  in  which  <  <  [h^  |  (see  for  example  eq. 
(127)  in  which  k  /  |k  |  &  J,  -.2  .  Hence,  k  in  the  exponent  of  eq.  (134)  can 
be  replaced  by  and  as  a  result  one  has, 
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where 
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The  substitution  of  eq.  (135)  into  eq.  (128)  (with  nt=2k^)  yields 
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or  upon  transforming  with  respect  to  t 
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Eqs.  (137)  and  (138)  still  constitute  rather  complicated  forms.  Severa  options 
are  open  for  their  further  reduction  of  which  two  are  discussed  below. 

1 .  First  Option 

The  averaging  process  over  V  can  be  carried  out  explicitly  un4er  the 
assumption  that 


Pk(£.  t  )  =  j  ds  V^(_r,  s) 


(139) 


is  a  multivariate  Gaussian  process.  The  Gaussian  character  of  cannot  be 

generally  justified  unless  t  is  large  compared  to  the  correlation  time  (T^)  of 

V  .  For  t  > >  t  p  is  Gaussian,  owing  to  the  central  limit  theorem,  regardless 
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of  the  precise  characterization  of  V  .  In  the  range  t  £  t  ,  p  is  Gaussian  only 
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re,  s'  =  si_s2  and  Cy(£,  s')  =  <  V^(r_,  s  j)  V  ( r_,  s^)  >  •  The  process  V_  has  been 
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presumed  to  be  transversely  homogeneous.  Consequently,  Cyis  a  function 
of  z  only.  If  we  also  assume  that  the  vertical  inhomogeneity  affects  pri¬ 
marily  the  variance  of  the  fluctuations  but  not  other  correlation  features 

such  as  t  ,  it  follows  that 
c 

Cy(r,S)  =  <V2(z)>  ty(s)  (141) 

where  fy(s)  denotes  the  correlation  coefficient  of  the  velocity  field.  The  sub¬ 
stitution  of  eqs.  (140)  and  (141)  into  (138)  yields, 
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or  alternatively. 


$E(uo)  =  M 


J dzL  J° dz2  y(Zj)  Y'  (z2)  $e(2ki’  z>  e^2kl  1-T1  Z  +  /  —  ~  Z 1  +  ^2> 


-h  -h 

2 

+  00  /* 

[ dT  e-jwT  -4k2  <V2(z2)>t 2  I  d|  iy(%) 

sin  2 

_|T 

J  o 

—  00 

_ 

143) 


where 


+  00 


1 

2tt 


/ 


,  A 

ds 


fy(3) 


.-A 


(144) 


represents  the  spectral-density  of  the  velocity  field.  If  we  further  assume 
2 

that  <V^  >  does  not  vary  (or  varies  slowly)  over  the  contributing  range  of  z, 
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one  observes  that  the  bracketed  terms  in  eqs.  (142)  or  (143)  are  independent 
of  z^  and  z^.  These  terms,  which  may  be  taken  out  of  the  respective  integrals, 
completely  determine  the  normalized  field's  spectrum: 


/s  A 

$e(ud)  =  M 


T 

oo  f 

j  dr  e~)m  e‘4kl2  <Vh  >  J  fV(s)  = 


A 

M 


-r00  r 

f  dTe-j“Ta-4kl2<V»>  T2/d5  {V<5) 

_ao  +® 

if 


v2 


(sin  2  \ 

4  ) 


with  M  selected  so  that 


duo  $E(uu)  =  1. 


(145) 


The  asymptotic  behavior  of  $E(u))  for  large  and  small  uu  (small  and  large 

T,  respectively)  can  be  readily  determined.  These  asymptotic  observations 

are  of  special  significance  since  they  are  independent  of  the  detailed  form  of 
A 

fv(s). 

(a)  The  limit  U)  o  (T>>TC)- 


The  major  contributions  to  the  ds  integral  in  the  exponential 
of  eq.  (145)  comes  from  the  range  §  ^  hence 


j  oo 

J  ds  (T  -  S)  fy(s)  w  |  T  I  J'  ds  fy(s)  =  I  T  I  Tc 


+  00 


/d 


jiDT  -4k,2  <v2  >  t  |t  |  = 


tE(°) 


0)  v2 


lim  $^(0))  =  M  /dx  e  JUM  e  H  ^  1  c  1 '  1  1+  (  ) 


U), 


0)“*O 


it  follows  that 


M  =  2k2  <V2  >  tc  $E(0)  =  1,  $E(0)  =  2/uj 


1 


(146) 


(147) 


(148) 
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and  the  defining  identities 


00 


ds" fy(s)  =  tt  *y(0) 


(JO.  =  4k-2  <v2  >  T 

1  1  H  C 


(149) 


(150) 


have  been  introduced. 

Comment ;  The  bound  character  of  the  vegetation's  scattering  centers 

(i.  e.  ,  the  finite  variance  of  their  displacement  in  the  limit  t  -*  ® )  implies 

a  DC  return  which  is  totally  ignored  in  the  above  sequence. 

(b)  The  limit  U)  -+  ®  (t  «t  ). 

c 

In  the  range  0  ^  s  <<  the  correlation  function  fy(s)  stays  virtually  constant: 

fy(s)  *  fy(0)  =  1  (151) 


hence, 


A 

lim  $E(tu)  = 
U)  -  00 


A 

M 


+® 
dr  e 


/ 


JOJT 


-2ki 
e  1 


<  V" 


2  ^ 

T  =  $E(0) 


Vrruu 


lTc 


U),  T 
1  C 


0) 


U), 


(152) 


The  transition  curve  connecting  the  asymptotic  segments  described  by 
eqs.  (147)  and  (152)  is  complicated  by  its  dependence  on  the  detailed  form 
of  fy(T).  f y( t  ) f  while  a  relatively  simple  characterization  of  the  coupling 
between  the  mechanical  motion  of  the  forest  scatterers  and  the  turbulent 
wind  is  an  unknown.  Little  reliable  data  exists  which  may  allow  its  determina¬ 
tion.  It  must  be  either  assumed  or  derived  on  the  basis  of  some  simple  (and 
probably  over  simplified)  models.  For  example,  one  may  regard  the  wind- 
forest  interaction  as  that  between  a  harmonic  oscillator  and  an  isotropic  and 
homogeneous  turbulent  wind  (M.  Labitt  and  R.  Yates,  private  communication). 
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Alternatively,  one  can  assume  that  regardless  of  the  driving  (wind)  forces 
the  resulting  motion  can  be  approximately  viewed  as  locally  harmonic  over 
the  relevant  time  interval.  This  option  is  discussed  in  the  second  part  of 
this  section.  The  form 


f  (t)  =  f  ,(t)  f  (t) 
v  d  o 


(153) 


where  f^(r)  an<^  fQ(T)  denote  decaying  and  periodic  functions  respectively,  is 
suggestive  (although  not  rigorously  defensible).  As  an  example,  selected  for 


its  relative  simplicity,  one  can  take  the  form 

It  I 


fy(T  )  =  e 


COS  n  T 


(154) 


which  permits  the  explicit  evaluation  of  the  integral: 

2 


/ 


(T  -  S)  fy(S)  it  = 


7^7" 

o 


-  cos  p+e  l/To  cos(p  +  n  It  |) 


(155) 


where 


cos  p  = 


-  2  2 

1-0  T 

i  +  oz°t2 

O 


o  <  (3  =  cos" 


i.a2T2 

T  <TT 


i+n2x 


(156) 


sin  p  = 


2Q  t 


“2  2 

l+cr  t 

o 


(157) 


Tc  is  defined  in  eq.  (149)  and  is  related  to  0  and  to  tq  via 


T  = 

c 


-2  2 
1  +  Q  T 

o 


(158) 


The  substitution  of  eq.  (155)  into  eq.  (145)  with  M  given  by  eq.  (148)  results  in: 
^Although  tc  as  given  by  eq.  (158)  should  not  be  strictly  interpreted  a,s  a  meas¬ 
ure  of  the  correlation  time  of  the  velocity  field  of  eq.  (154),  it  is  convenient 
to  preserve  the  definition  given  by  eq.  (149). 
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'♦E(“)  =  2kl2  *E<0)<  v2>  Tc 


+  00 

/ dT  e 


e-4kf<V2>  to  Tc 


-  cos  3+  e 


“  T 


/  T 


O  COS  (3+0  |  T  I  ) 


(159) 


One  immediate  conclusion  stemming  from  eq.  (159)  is  that  even  for  the  simple 
choice  (154)  no  closed  form  expression  for  1^(uu)  is  available.  The  Fourier 
transform  (159)  was  evaluated  on  a  computer  and  examples  of  these  results 
are  shown  in  Fig.  10.  The  theoretical  result  given  by  the  solid  curve 
(Qtq  =  5,  ut)^To  =  0.5)isvery  similar  to  measured  UHF  clutter  spectra  with 
moderate  winds. 

2.  An  Alternative  Option 

The  forthcoming  analysis  starts  once  again  with  eq.  (138).  However, 

instead  of  presuming  the  process  V  (r_,  t)  to  be  normal  we  now  postulate  a 

H 

locally  harmonic  temporal  dependence: 


V^(r_ ,  t)  =  vm(rJ  cos  [Q  (£)t+cp(ij  ] 


(160) 


where  V  (r),  Q(r_)  and  cp(r_)  are  random  functions  of  r_  but  are  time  independent. 
The  relation  (160)  is  not  expected  to  constitute  a  valid  description  over  pro¬ 
longed  periods.  The  resulting  spectral  predictions ,  therefore,  are  expected 
to  be  inaccurate  near  DC.  In  this  respect,  the  model  to  be  pursued  in  this 
section  can  be  regarded  as  complementary  to  that  leading  to  eq.  (147).  It 
follows  trivially  that, 

T 

v  (£,  s)  ds  =  p  (r)  [sin  n  T  +  cp)  -  sin  cp  ]  ( 1 6  1 ) 

h  m  — 

o 


i 

f 


where  V  and  p  =  V  /Q  represent  local  maxima  of  V  and  the  displacement, 
m  rm  m  H 


Vanishing  of  the  field  as  t -*  00  implies  unbound  scatterers.  We  are  aware 
that  the  scatterers  are  indeed  bounded  but  here  have  chosen  to  disregard  this 
feature  since  the  loss  of  a  measure  of  the  stationary  backscatter  (impulse  at 
UU  =  o)  is  not  a  major  penalty. 
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Fig.  10.  Clutter  spectral  density. 
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respectively,  If  we  assume  p  ,  0  and  cp  to  be  statistically  independent  (not 
entirely  unreasonable  since  Q  and  cp  are  expected  to  be  essentially  independ¬ 
ent  of  wind  conditions,  while  a  strong  wind  dependence  is  anticipated  for 

n  ),  one  obtains 
Km 


IE(u>) 


+  00 

*</ dT  e 


-jU)T  +  j2k 


./ 


ds 


Vr- 


s)  y 


<e 


+  eo 

•11  ,,  i  \  ■  I  j  -juu  t  +  j2ki  pTY,(r)  sin  (Q  t  +  cp)  v 

-j2kL  pm(_r)  sin  cp  /  dT  e  J  J  1 

J  o 


(162) 


m 


cp 

Q 


where  uu )  is  defined  by  the  expression  to  its  right,  and<^ 

cpm 

n 

indicates  sequential  averaging  over  the  distributions  of  p  ,  cp,  and  Q  (the 

actual  averaging  sequence  is  arbitrary).  The  substitution  of  the  expansion 

+  00 


ej2kiPmsin(QT+cp)  =  £  ^  (2k  pm)  ejn(nT+cp) 

n=  -oo 

into  eq.  ( 1 6 2 )  and  carrying  out  the  t-  integration  results  in 

00 

y<»>  =  2n  <e-i2kl  pm  =  “  v  £  Jn  (2k;  pm)  <>*  6<u>-  »fll  > 


(163) 


(164) 


n=  -oo 


m 


cp 

n 


The  averaging  over  cp  is  readily  performed  if  we  presume  a  uniform  dis¬ 
tribution  within  the  interval  -  tt  <  cp <  TT  : 

TT 

<e-j2kl  Pm  Sin  "P  +  JncP>  =  J-  f  dcpe-j2klPmsinCP+jncP^J  (2k,p  ) 

cp  2tt  J  n  1  m 


(165) 
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Eq.  (164)  is  reduced  to 

”  2 

Ir(uo)  =  2tt  Y  <J  (2k,  p  )\  /'.  n.s 

E'  7  ^  \  n  v  1  m;  /p  \0(uu-  = 

n=  -® 

00 

2tt  <Jq2  (2kL  pm)>  6(uo)  +  2tt  y  <JIf(2k1Pm)>  <  6  (u)-nQ  )  +  6(uu+nn>n  ( 166) 

pm  n=  1  pm 

Let  P(fl)  denote  the  density  function  specifying  the  Q  distribution  (its  general 
properties  are  discussed  below).  The  0  averaging  leads  to 

IE(U>)  =  2tt  <Jq2  (2k1Pm)  >  6  (U))  + 

00 

2tT^.  n  <Jn(2kipm)>p  [P  ( )  +  P  (- -!£-  )  ]  (167) 

n=  1  Km  L  '  n  v  n  J 

It  can  be  immediately  observed  that  I^(ud)  is  symmetric  about  uu=  o,  regardless 

of  the  precise  nature  of  P(Q). 

The  task  of  averaging  over  is  more  difficult.  However,  before  turn¬ 
ing  to  it  we  discuss  some  of  the  properties  of  eq.  ( 1 6 7)  which  are  relatively 

2 

insensitive  to  the  term  <J  (2k,  p  )>  n 

n  v  1  Km  Pm 

P(Q)  is  not  known.  However,  since  it  essentially  describes  the  distribu¬ 
tion  of  resonance  frequencies  (positive  by  definition)  of  the  ensemble  of  forest 
constituents  it  must  vanish  in  the  limits  fi  -*  »  and  0  £  o.  It  must  contain 
at  least  one  maximal  value;  for  simplicity  let  us  assume  that  P(Q)  possesses 
a  single  maximum  (say  at  fl  =  Q^). 

The  physical  interpretation  of  eq.  (167)  is  straightforward.  It  contains 
two  distinct  parts.  The  first  (in  our  case  proportional  to  6  ( uo ) )  represents 


57 


that  part  of  the  received  power  which  did  not  undergo  spectral  modifications. 

Its  spectral  content  is  identical  to  that  of  the  transmitted  signal,  although 

2 

it  experiences  a  power  reduction  (represented  by  the  coefficient  <JQ(2k^  Pm)> ) 
as  the  power  conversion  process  into  other  spectral  constituents  (represented 
by  the  second  term  in  eq.  ( 1 6 7)) becomes  more  efficient.  The  efficiency  of 
this  conversion  process  depends  primarily  on  the  (random)  parameter  Cy¬ 
linder  no  wind  conditions  (k.  p  =  0  )  all  the  power  is  contained  in  the  first 
(DC)  term.  The  second  term  vanishes.  As  k^  pm  increases  (with  the  rising 
magnitude  of  wind  velocity),  the  DC  term  decreases  in  magnitude,  the  second 
term  increases,  consistently  with  the  anticipated  conservation  of  energy. 

The  singular  character,  of  the  DC  term,  stems  directly  from  the  presumed 
monochromatic  nature  of  the  transmitted  signal.  Eq.  (167)  stays  valid  if 
the  actual,  distributed  spectrum  of  the  transmitted  signal  replaces  the  6  - 
function,  provided  that  it  is  sharply  peaked  (on  the  scale  of  0q). 

The  spectral-density  $^(uu)  can  be  obtained  now  by  the  substitution  of 
eq.  (167)  into  eq.  (138).  However,  because  of  the  complexity  of  the  result 
we  presently  assume  p^  to  be  statistically  homogeneous  in  all  directions. 
Consequently,  I^(tl))  is  independent  of  z^  and  z^  and  represents  by  itself  the 
desired  spectral  density. 

2 

We  now  turn  to  the  estimate  of  the  coefficients  <J  (2k,  p  )>  in 

n  1  lrm  p 

rm 

the  limits  of  weak  (k,  p  <<  1)  and  strong  (k,  p  >  >  1)  winds. 

(a)  The  limiting  case  k^  Pm  <<1. 

The  Bessel  functions  may  be  replaced  by  the  term  in  the  corre¬ 
sponding  power  series  expansion: 
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2  <^klPrr,)2n>D 

<J  2k.  p  >  tv  - 5 - — 

n  '  l^m'  p  ,  ..2 

rm  (n!) 

| 

(168) 

The  series  is  a  rapidly  convergent  one  and  for 

<(k.  p  )2  >  >>  4  <(k.  p  )4> 

v  1  m'  p  8  1  m'  p 

m  rm 

(169) 

it  is  properly  represented  by  its  first  term: 

Ie(«m)  «  lEo)(t«)  +  2n<(k1  pm)2  >  [Pu>)  +  P(-uo)] 

pm 

(170) 

where  lj*°^  represents  I  _( uu )  measured  under  no  wind  conditions  (p 

iii  A  iii  m 

In  eq.  ( 167),  (uu)  is  identified  as  2tt  6  (uu) . 

iii 

f  o). 

The  presumption  of  validity  of  the  above  model  indicates  the  feasibility 

of  the  experimental  determination  of  P(uu)  from  the  measurement  of  I 

:E(w) 

under  weak  wind  conditions.  From  eq.  (170)  one  has 

P(ud)  «C  [Ie(uj)  -  I^o)  (UJ) ]  u  (id) 

00 

(171) 

where  C  is  selected  consistently  with  the  requirement:  C  J [i^uu)  “ 

>)(uu)]duu=l. 

o 

The  experimental  determination  of  P(uu)  via  eq.  (171)  is  significant  owing  to 

the  anticipated  insensitivity  of  P(uu)  to  varying  wind  conditions.  The  c 

ieter- 

mined  P(uu)  may  be  utilized  in  eq.  (167)  for  relatively  arbitrary  winds 

• 

.  (b)  The  limiting  case  k^Pm  >>1.  An  exact  evaluation  of 

<  J2  (2k,  p  )> 
n  Pm'  p 

m 

to 

It  will  be  found  advantageous  to  represent  the  Bessel  functions  in  eq. 

(167) 

by  their  integral  representation.  We  have, 
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<  J  2  (2k.  p  )>  h  i  <  [(H(1)  (2k.p  )  +  H(2)  (2k,  p  )]2> 
n  1  rm  p  4  L  n  v  1  Krrv  n  v  1  i 

m 


m 


2  2 

i<H(1)  (2k  p  )  +  H(2)  (2k.  p  )+2H(1)(2k.  p  )  H(2)  (2k,  p  )>  (172) 

4  n  1  m  n  v  1  Km  n  v  1  Km  n  v  1  m  p 


m 


where  the  Hankel  function  is  defined  via 


(2k  p  )  =  -  f  e_j2kl 

n  1  Km'  rr  / 


pm  sin  cp  +  jncp  dcp 


(173) 


1 

2 


m 


and  the  paths  are  illustrated  in  Fig.  11.  For  simplicity  we  assume  p 

2  _  _  2  1/2 
to  be  normally  distributed  with  a  mean  p  and  a  variance  o  =  <(p  -p  )•  > 

7  Km  p  m  m 

If  we  further  define  the  parameters; 

1/2 

P 

we  find, 


L  =  [  (2kx  p'm)2  +  (2k2  o2  )2]  ,  sin  d  = 


2k.  p 

1  rn  .  tt 

- -  ,  o  <  Of  <—  (174) 


<2)2 

< H  (2k.  p  )  > 
n  '  1 rm  p 


m  tt 


Jj.  J  dcpj  /dtp2  ejn(Pl+<P2>  e-LQItP!,  »2)  (175) 


1 

2 


and 


<H(1)  (2k.  p  )  H(Z)  (2k.  p  )  > 
n  v  1  Km  n  v  1  rm  p 


m 


=  fa  fa  .*<»,♦,*) -U*,,.,*) 


C1  C2 


(176) 


where 

2 

CHcp^ ,  cp^ )  =  j  sln  O'  (sin  cf>j  +  sin  cp^)  +  cos  a  (sin  cp^  +  sin  cp^)  (177) 
Eqs.  (175)  and  (176)  are  now  evaluated  asymptotically  (with  L  taken  to  be 
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Fig.  11.  The  path  defining  Hn(-|-). 
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the  large  parameter)  via  the  method  of  steepest  descent.  The  saddle  points 


cp^s  and  defined  by  the  equations 


tois.cp2s) 


=  0  = 


2  s 
TT 


n 


a  cpls  “  d  cp 

are  readily  found  to  be  situated  at  either  +~(  on  )  or  at  -  —  (along 
Eq.  (175)  yields  the  saddle-point  contribution. 


(178) 


(V 

.  u '2 '  /0,  >  .  2n 

<  H  (2k  P  >  ^  -r—  _ - rr 

n  1  m  p  L  f  -v  2 

m 


jn(cp  +  cp  ) 

is  2s  -LQ(cpls,  cp2s)  +  Q  (  J_} 


d  Q  d  Q 


b  cp2  a  cp2 


172 

V^s 

*2=,%s 


Ztt 


TV 


jnTT  -  L  (4  cos  <y  +  2j  sin  a  ) 
±  j  sin  o'  -  4  cos  a 


+  O  (-V) 


where  the  (±  )  corresponds  to  (^) •  Similarly, 


(179) 


<H(1)(2k.p  )  H(2)  (2k.  p  )>  ~  —  +  O  (  — k 

n  '  1  Km'  n  v  1  Km  p  L  sin  q-  T  2 

m  L 


(180) 


Hence,  from  eq.  (172) 


<  J  (2k.  p  )  > 
n  v  1  lrn/  p 


m 


-  f-  -■  ■ 

L  I  sm  a 


1 

r  .  2  ~7  2  T72 

[  sm  a+  16  cos  a  \ 


e  COS  ^  cos  (2L  sin  a  -  nrr  +  0)1,  |3  =  tan  ^  [  -j  tan  <y  ]  (181) 


A  substantial  simplification  occurs  if 

4L  cos  a  =  8kf  a ^  >>  1 
1  P 


(182) 
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The  last  term  may  simply  be  omitted  owing  to  the  exponential  decay  resulting 
in 


<  J  (2k.  p  )> 
n  1  rm  p 


m 


TT 

2k.  (5 
1  m 


(183) 


which  is  independent  of  n.  Hence,  in  the  limit  cjp>>  1  the  normalized 
spectrum  ( 1 6 7)  now  reduces  to 


2  oo 

!E<“>=  Tl-  !6(u,)+  ^  i  +  1  (  (l84» 

K1  pm  (  n=l  ) 

The  asymptotic  expansions  (179)  (180)  presumed  not  only  the  largeness  of  L 

with  respect  to  unity  but  also  with  respect  to  n.  In  the  case  where  the  series 

( 1 6 7 )  does  not  converge  at  a  sufficiently  rapid  rate  the  results  (179)  (1$0)  must 

be  reconsidered.  The  asymptotic  analysis  must  allow  for  large  values  of  n. 

The  saddle-point  condition  and  the  corresponding  contributions  are  modified 

as  indicated  next.  Let, 


6  =  jn  (cp^  +  cp^)  -  L  [j  sin  or  (sin  cp^  +  sin  cp^)  +  cos  afsincpj  +  sincp^)2] 


(185) 


hence. 


d  Q 


d  cp 


jn  -  L  [  j  sin  a  cos  cp^  +  2  cos  a  (sin  cp^  +  sin  cp^)  cos  cp^  ] 

2  2 


(186) 


and 


2  A 

a  q 

d  Cfj 


2  2 

L  [  -j  sin  a  sin  cp^  4-  2  cos  a  (cos  cp ^  -  sin  cp^  -  sin  cp^  sin  cp^ )  ]  (187) 


The  saddle-point  conditions  are  defined  by 

a q  |  =  a6  | 

’is.  *28  '  **2  Ss,  ^2s 


(188) 
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We  specifically  look  for  solutions  of  eqs.  (188)  subject  to  the  constraint 

cpis  =  -cp2s  (189) 

which  selects  the  unique  combination  of  cp  and  cp?  giving  rise  to  contributions 

J.  s  L*  s 

which  are  not  exponentially  small.  Such  a  solution  indeed  exists  and  is  given 
(via  eq.  (186))  by 


jn  -  jL  sin  a  cos  cp  =  0 

Is 


(190) 


or 


cos 


n 


Kls  L  sin  a  1 
2 


sin  cp,  =  ± 

Ms 

2 


1  - 


(-r - " -  )2  (191) 

L  sin  a 


The  substitution  of  eq.  (191)  into  (187)  results  in, 
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and 


.  2  2 

a  q  a  q 

N  2  7  2" 

O  cp  1  d  cp2 


.,2  -  2  2,4 

4k,  p  -n  +n 
1  m 


=  L 

'Pis  =  -fz. 
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Therefore, 


<J  (2k.  p  )> 
n  '  1  rm  p 
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2  -  2  2  4  /  \4 

p  -n  +n 
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<f) 
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(194) 
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and  substitution  into  eq.  ( 1 6 7 )  results  in 


IE(U)) 


TT 


kl  P™ 
1  m 


6  (ud)  +  2tt  £ 


P(  ®)  +  P(Jfi) 

n  n 


n=l  2  -  2  2  4  /  ap  V 

nV4kl  pm  _n  +n  (ipj 


195) 


An  exact  evaluation  of  <  J  (2k.  p  )>  is  possible  for  a  presumed 

n  1  p  r  ^ 

m 


Rayleigh  distributed  variable  p  .  One  has 


<J  (2k  p  )> 
n  v  1  Krrv  p 


m 


00 

/ 


dp  P(p  )  J  (2k,  p  ) 
m  m  n  1  Km 


—  f  dp  p  J2  (2k  p  )  e"  A  2  =  I  (4k2  £  2)  e“4kl  Pm  (196) 
,  ^  /  Km  m  n  1  2  p  n  1  m 

^  m 

m  o 


A. 

P 


where  p  denotes  the  value  of  p  at  the  peak  of  P(n  ). 
m  m  r 

Hence,  via  eq.  (167) 

2  2  oo  £  ^  2 

IE(»)  =  2n  Io(4k12Sm2e-4kl  6(0,1  +  £  i  !„(  4k,2  gjl  e‘4kl  °m[p,  1?)  +  P(.  |)  ] 

(197) 

which  reduces  to  the  results  given  in  eqs.  (170)  and  (184)  in  the  respective 
asymptotic  limits. 

The  clutter  spectrum  given  by  eq.  (197)  again  indicates  a  DC  component 

at  uo  =  o  plus  continuous  AC  spectrum  with  the  appropriate  interchange  between 

2  2 

these  components  to  conserve  energy  as  the  factor  4k^  is  varied.  I^(u)) 

as  given  by  eq.  (197)  is  normalized  for  unit  total  clutter  power.  Hefnce, 

■foo 


p  =  p  +  p 
T  DC  AC 


■if 


IE(u' )  dll)  =  1 


(198) 


The  DC  clutter  power  is 
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(199) 


„  T  2  A  2,  -4k.  ft  2 

Pnr  =  I  (4k,  ft  )  e  I'm 
DC  o'  1  Km 

Hence,  the  AC  clutter  power  is  given  by 

PAC  =  PT  -  PDC  =  ‘  -  :o  <4kl2  C  »  ' 


(200) 


Note  that  this  is  true,  independent  of  the  specific  form  of  the  density  function 
of  self -resonant  frequencies  P(fi).  Figure  12  is  a  plot  of  the  ratio  of  AC  clutter 
power  to  DC  clutter  power  as  a  function  of  the  parameter  ^  /\ .  For  large 


m 


displacements,  this  ratio  is  proportional  to  p^/X. 

Figure  13  depicts  a  measured  clutter  spectrum  at  UHF  (X  =  0.69  m)  under 
low  wind  conditions.  All  of  the  measured  spectra  were  characterized  by  a 
resonant  peak  at  approximately  0.4  Hz.  The  linear  slope  (db  vs.  log  fre¬ 
quency)  of  the  observed  low  wind  spectra  above  this  peak  suggests  a  power 
law  representation  for  P(fi)  of  the  form 

a  (n  /Q  )P 


p(n) 


(201) 


1  +  b  (0/Qo)q 

(see  eq.  (171))  .  A  good  fit  to  the  tail  of  the  low  wind  spectrum  was  obtained 
with  parameters  chosen  such  that  q  -  p  ^  5.7.  Alternative  choices  for  P(Q), 
such  as  a  Rayleigh  distribution  or  the  Maxwell  distribution  included  in  Fig. 

10,  can  not  accurately  produce  the  observed  low  wind  clutter  spectra. 

The  power  law  distribution  which  matched  the  low  wind  UHF  spectra  was 
then  used  in  eq.  (197)  to  predict  clutter  spectra  at  both  UHF  and  L.-band 
carrier  frequencies  for  higher  wind  conditions.  These  results  are  illustrated 
by  Fig.  14.  The  experimental  data  at  the  two  frequencies  was  obtained  simul¬ 
taneously  from  a  common  range  cell.  It  is  seen  that  the  power  law  distribution 
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Fig.  12.  AC-to-DC  clutter  ratio  for  Rayleigh  distributed  displacements. 
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Fig.  13.  Low  wind  clutter  spectra.  Deflections  Rayleigh  distributed. 
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Fig.  14.  Moderate  wind  clutter  spectra.  Deflections  Rayleigh  distributed. 


t 
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with  a  deflection  parameter  of  =  3.7  cm  which  quite  accurately  matches  the 
observed  L-band  spectrum  also  predicts  reasonably  well  the  observed  UHF 
spectrum. 

The  Maxwell  distribution  of  P(Q),  which  was  seen  to  be  a  poor  choice  under 
low  wind  conditions,  yields  much  more  reasonable  results  with  stronger  winds. 
This  suggests  that  the  higher  wind  spectra  are  less  sensitive  to  the  distribution 
of  self -resonant  frequencies  of  the  scatterers. 

It  should  be  observed  that  the  distribution  of  characterized  by  the  para¬ 
meter^  may  depend  upon  the  frequency  of  the  incident  radiation.  As  the  fre¬ 
quency  increases ,  smaller -scale  scattering  centers  start  to  effectively  partici¬ 
pate  in  the  scattering  process.  It  is  anticipated  that  these  smaller  constituents 
will  move  with  larger  amplitudes  and  thus  thend  to  increase  .  The  trend  is 
reversed  as  the  frequency  decreases.  While  the  qualitative  observation  is 
simple,  no  satisfactory  quantitative  procedure  has  been  established. 
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APPENDIX 


The  objective  of  this  Appendix  is  to  prove  the  validity  of  the  relation: 

Var  P=<[P-<P>]2=<P>2  (Al) 

r  L  r  r  r 

under  the  constraints  of  single-scatter  theory  and  the  presumed  largeness 
(compared  to  the  correlation  length)  of  the  scattering  volume. 

Let 

E  1  =  Re  [  A  Eg.  ] ,  E2  =  Im  [  A  Egi  ]  (A2) 

The  largeness  of  AV\  in  eq.  (118)  assures  (via  the  central  limit  theorem)  the 
normal  character  of  A_Eg^  (which  is  therefore  completely  defined  by  its  two- 
point  correlations).  Upon  forming  the  expression  <4  E  .  ■  A  E  .  >  (via  eq. 
(118)),  transforming  the  resulting  double  integral  (say  over  r^  and  r^)  into 
the  center  of  mass  (R  =  (r_  +  v_^))  anc*  relative  (L  =  L\  “  JL?)  coor<^lnate 

systems  and  further  separating  the  horizontal  coordinates  from  the  vertical, 
one  obtains: 
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It  is  the  bracketed  term  in  eq.  (A3)  which  is  of  special  interest,  since 
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Hence,  the  largeness  of  A  is  measured  on  the  scale  of  the  wavelength  of  the 
incident  radiation  (rather  than  the  correlation  distance  characterizing  e  (£,  t). 
It  follows  that, 

<A_Eg^  (£,  t)  -  A_Eg^  (r^,  t+T  )>  =  0  (A5) 

or 

<_Ej  (£,  t)  •  JEj  (r/,  t+  t )  >  =  <  Ez  (r,  t)  •  (r_‘,  t+T)>  (A6) 

and 

<Ej  (£,  t)  •  _E2  (£r,  t+T)>=  -<  E1(£/,  t+T)  •  _E2  (£-  t)>  (A7) 

Although  one  can  continue  with  the  treatment  of  the  two-point  statistics, 
we  restrict  ourselves  to  the  discussion  of  the  one  point  statistics  relevant  to 
the  proof  of  eq.  (Al).  With  r_  =  r  '  and  t  =  0  in  eqs.  (6,  7)  it  follows  that: 


<E12>  =  <E22>  =  0  E  ’  <E1  '  E2>=  0 


Hence,  the  joint  distribution  is  given  by 
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Or  in  terms  of  the  amplitude  and  phase,  with 


A  =  [Ex2  +  E2  ]1/2,  cp  =  tan-1 

one  obtains 
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For  the  Rayleigh  distribution  (All)  it  follows: 

2  2  4  4 

<  A  >  =  2  Og  ,  <  A  >  =  8  Cg 
Hence  (with  S  denoting  the  Poynting  vector), 

Var  S  =  <[S-<S>]2>  =  ±  —  <  [  A2  -  <  A2>  ]  2>  = 


4 

where  eq.  (Al  2) 


* 


^  [<a4>  -  <a2>2  ]  =—  c  4  =<s>; 

Li  E 
p  o 


was  utilized.  Eq.  (Al)  results. 
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